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Abstract. We present a quantitative estimate of the accuracy of the calculated diffusion coefficients, by comparing 
predictions of solar models with observational data provided by helioseismology. By taking into account the major 
uncertainties in building solar models we conclude that helioseismology confirms the diffusion efficiency adopted in 
SSM calculations, to the 10% level. 
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1. Introduction 

In recent years inclusion of elemental diffusion has been an essential ingredient of stellar evolutionary codes for achieving 
agreement between predicted and helioseismic values of properties of the convective envelope, see e.g. Cox et al.(1989), 
Bahcall & Pinsonneault (1995), Bahcall et al. (1997), Ciacio et al. (1996), Degl’Innocenti et al. (1997), Richard et 
al. (1996), Turck-Chieze et al. (1998), and Fig. 1. The success of solar models with diffusion, and the corresponding 
failures of models that neglect diffusion, suggest that the diffusion process has been properly treated. The goal of 
this paper is to present a quantitative estimate of the accuracy of the calculated diffusion coefficients, by comparing 
predictions of solar models with observational data provided by helioseismology. 

Until recently, few stellar evolution codes included diffusion. Apart from practical difficulties, the omission of 
diffusion was justified since elemental diffusion occurs on a large time scale (larger than 10^^ yr. to diffuse a solar 
radius under solar condition), which should imply that the effect on stellar structure are small. Nevertheless, the 
precision required for calculating solar neutrino fluxes and for matching the accuracy of helioseismic observations is 
so great that diffusion has to be taken into account, see e.g. Bahcall & Loeb (1990), Proffit (1994). 

Most importantly, diffusion together with gravitational settling hide below the convective envelope a significant 
fraction (about 10%) of the initial helium and heavy metals, and this effect is crucial when comparing solar models 
with the observed properties of the convective envelope (e.g.the present heavy element abundance and helium fraction, 
this latter inferred by means of helioseismology). 

One of the most detailed calculations of diffusion coefficients has been presented by Thoul et al. (1994), hereafter 
TBL. It is based upon exact numerical solutions of the fundamental equations for element diffusion and heat transfer 
discussed in Burgers (1969). Good agreement exists between the diffusion rates computed by Thoul et al. and the 
results obtained with a very different treatment of Michaud & Proffit (1993). The accuracy of the TBL coefficients is 
estimated to be about 15% by comparison with the results of other authors, see Bhacall and Pinsonneault (1995) and 
refs, therein. 

However, some warning is needed: the separation of elements by diffusion might be inhibited by the presence of 
mixing in the solar interior, see Schatzman (1969), Lebreton & Maeder (1987), Pinsonneault et al. (1989), Pinsonneault 
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et al. (1990), Chaboyer et al. (1992), Chaboyer et al. (1995). We also remind the anomalous lithium depletion of the 
sun with respect to meteorites, signalling some deficiency of present standard solar models. 

Clearly quantitative observational constraints on the calculated diffusion coefficients would be welcome. The pur¬ 
pose of this paper is to discuss the constraints posed by helioseismology. 

As previously remarked recent Standard Solar Models (SSMs) which include diffusion, as well as accurate equations 
of state and updated opacity tables, are in good agreement with helioseismology. In essence we address the following 
question: which variation of the diffusion efficiency is allowed without spoling the agreement with helioseismic data? 

In the next section we study how the predictions of solar models are affected when the diffusion efficiency is changed 
and we summarize the relevant helioseismic information. 

In sect. 3 we discuss the comparison of helioseismic data with the predictions of solar models, when the diffusion 
efficiency is varied. 
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Fig. 1. The photospheric helium mass fraction Yph and the depth of the convective zone {Rb/Ro). The dotted rectangle is 
the Icr interval allowed by helioseismology from Degl’Innocenti et al. (1997); the open circles are the predicted values by solar 
models without diffusion, from top to bottom correspond to Dar & Shaviv (1996), Christensen-Dalsgaard et al. (1993), Richard 
et al. (1996), Turck-Chieze & Lopes (1993), ProfBt (1994), Bahacall & Pinsonneault (1995), Ciacio et al. (1996). The full squares 
are the predicted values by solar model with helium diffusion, from top to bottom correspond to Christensen-Dalsgaard et al. 
(1993), Bahacall & Pinsonneault (1992), Proffit (1994). The full circles are the predicted values by solar model with helium and 
heavy elements diffusion, from top to bottom correspond to Richard et al. (1996), Cox et al. (1989), Proffit (1994), Bahacall & 
Pinsonneault (1995), Christensen-Dalsgaard et al. (1996), Turck-Chieze et al. (1998), Bahacall et al. (1998), Ciacio et al. (1996), 
Dar & Shaviv (1996). 



2. The diffusion efficiency and solar models predictions 

For our standard solar model calculations we use the diffusion coefficients calculated following TBL. For each element 
one has to consider three different coefficients {Ap, Ap, and Ap following the notation of TBL), each depending on the 
value of several chemical and physical variables (p, T, X, Zi...) in the solar interior, so that in principle one has really 
many parameters which could be varied for numerical studies. All the diffusion coefficients, however, are approximately 
constant along the solar profile apart from the innermost core, see fig. II in TBL. 
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For simplicity, we shall introduce just one diffusion efficiency parameter d and we shall study solar models obtained 
by rescaling the TBL coefficients by the factor d, assumed to be the same for all elements at any point in the solar 
interior. 

By numerical experiments with our evolutionary code FRANEC (Cieffi & Straniero (1988)) we have determined 
solar models (i.e. stellar structures with luminosity Lq = 3.844-lO^^erg/s, radius Rq = 6.9599-10^°cm and photospheric 
abundance ZfXph = 0.0247 at age = 4.57 Gy for a mass Mq = 1.989 - 10^^ g) obtained for different values of d. 

By means of helioseismology one can reconstruct accurately three independent properties of the convective envelope: 
its depth Rb, its density at the bottom pb and the photospheric helium abundance Yph- The resulting values together 
with estimated Icr uncertainties are shown in Table 1 together with SSM predictions from the model with helium and 
heavy element diffusion of Bahcall and Pinsonneault (1995) hereafter BP95. 


Table 1. The solar observables Qi used in our analysis, their predicted values QssMi from BP95, the helioseismic determination 
Qobi with their la error AQi, from Degl’Innocenti et al. (1997) and the calculated coefficient a, corresponding to Eq. 1. 


Qi 

QsSMi 

Qei 

AQ, 

OtQ 

Yph 

0.247 

0.249 

0.0035 

-0.091 

Rb/ Rq 

0.712 

0.711 

0.0014 

-0.016 

pb [g/cm’^j 

0.187 

0.192 

0.0018 

0.14 


By calculating these observables for our solar models, we have studied their dependence on d. For each observables 
Q we introduce a parametrization of the form: 

Q(d) = QssMd°‘^ (1) 

The coefficients a are collected in Table 1. 

The numerical results can be qualitatively explained by simple considerations. Consider as an example the case of 
increased diffusion efficiency (d > 1), with respect to the SSM (d = 1): 

i) As a larger fraction of helium is hidden below the convective envelope, whereas the initial helium abundance is 
essentially fixed by the present luminosity, this implies a lower photospheric helium than in the SSM; 

ii) the metal abundance below the convective envelope (for a fixed photospheric abundance) is increased so that opacity 
increases and convection starts deeper and consequently at higher density. In addition from helioseismology one can 
reconstruct the sound speed profile of the solar interion with a la accuracy of about 0.15% in the intermediate region 
(0.2 < R/Rq < 0.6). 

This piece of information, however, is not particularly illuminating for studying of diffusion. As already remarked 
in Degl’Innocenti et al. (1997), solar models without diffusion can have sound speed profile quite consistent with the 
helioseismic one. 


3. Diffusion efficiency and the properties of the convective envelope 

We now determine the acceptable range of d by requiring that the three independent observables of the convective 
envelope, Rb, Pb and Yph-, are predicted within their helioseismic ranges, by using Eqs. (1) to determine the dependence 
of these properties on d. 

We remind that there are at least three major uncertainties in building standard solar models that also have the 
potentiality of affecting the three properties under investigation, and, therefore, that could interfere with/hinder the 
effect of d'. the astrophysical factor Spp, the solar opacity k, the heavy element abundance ( = ZjX. We shall add all 
these effects one after the other, and determine a range of d that takes into account these uncertainties. 

We start by defining a as: 


‘w = E 

i=l,3 


Qiidf) Qqi 

Ag) 


( 2 ) 


where the sum include the three independent observables of the convective envelope, Qi{d) are computed by using 
Eqs. (1), and Qqi are the helioseismic values reported in Table I together with the Icr errors XQi. The value ^^{d = 1) 
indicates how well the SSM reproduces these helioseismic properties. The first row of Table 2 shows the good agreement 
between BP95 and helioseismology (x^/dof = 8.61/3). 
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If we use d as free parameter (second row of Table 2), we find the following best fit value (x^/dof = 2.92/2) and 
1(7 range: 

d= (1.14 ±0.06), (3) 

i.e. helioseismology suggests a diffusion efficiency slightly larger than the SSM estimate. 


Table 2. Deviations from standard diffusion allowed by helioseismic measurements. The first five columns show whether the 
parameter is kept fixed (F) at its SSM value or it is allowed to vary (V) as a free parameter. The sixth column shows the 
resulting per degree of freedom. The last two columns show the best fit value for d and its Icr error. 


K 

C 

Spp 

d 

xVdof 

C^Best 

Ad 

F 

F 

F 

F 

8.61 /3 



F 

F 

F 

V 

2.92 /2 

1.14 

0.06 

F 

F 

V 

V 

1.81 /2 

1.07 

0.09 

F 

V 

V 

V 

0.98 /2 

1.05 

0.09 

V 

V 

V 

V 

0.75 /2 

1.08 

0.11 


3.0.1. Uncertainties on Spp 

A conservative estimate of the uncertainty is provided by the range of the published results (NACRE coll. (1998)), 
whereas a Icr estimate has been provided in Kamionkowski & Bahcall (1994); we shall use AS'pp/S’®®'^ = 0.05/3 at Icr 
(5% is the “3cr error” estimate). The dependence of Qi on Spp has been determined numerically in Degl’Innocenti et 
al. (1998). By redefining a suitable x^{d,Spp): 


\d,Spp) = 


Qiid^Spp^ Qqi 

Aft 


Spp *S^pp,SSM 


AA 


pp 


( 4 ) 


we find now: 


d= (1.07 ±0.09), 


( 5 ) 


3.0.2. Uncertainties on k and ( 

The heavy element abundance C and the solar opacity k are known with a conservative accuracy of about 10%, see 
Bahcall & Pinsonneault (1992), V. Castellani et al. (1997). Therefore, our Icr relative error estimate will be 0.1/3. The 
dependenceWe remark that we are considering a constant rescaling of opacity along the solar profile, of Qi on k and 
C has been determined numerically in DeglTnnocenti et al. (1998). In this case, the relevant is: 


c A f Qiidj Spp,Q k) Qoi\ 

X {d,Spp,Q,K) - I 

i ^ ' 

/ Spp - S'pp.SSM \ ^ 

V A5pp ; V ) 


( 6 ) 


We find a small change of the best fit value and, of course, an increase of the Icr range of d\ 


ci= (1.08 ±0.11). 


( 7 ) 
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3.0.3. Solar model “theoretical uncertainties” 

At last we try to estimate how much our results could depend on having used BP95 as reference standard model. To 
this end, we consider one of the standard solar models (models that include all the state-of-the-art solar physics), whose 
helioseismic properties differ the most from BP95 and, consequently, fit less well the experimental data. We repeated 
the above-described analysis by using FR97 (Ciacio et al. (1996)) as standard solar model. When all parameters are 
varied, the Icr range and best fit value, c/. Eq. (7), become: 

d= (1.13 ±0.12), (8) 

in substantial agreement with Eq. (7). 

In conclusione, helioseismology confirms the diffusion efficiency adopted in SSM calculations, to the 10% level (Icr 

C.L.). 
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